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Abstract 

The ferromagnetic Ising spins are modeled on a recursive lattice constructed 
from random-angled rhombus units with stochastic configurations, to study the 
magnetic properties of the bulk Fe-based metallic glass. The integration of spins 
on the structural glass model well represents the magnetic moments in the glassy 
metal. The model is exactly solved by the recursive calculation technique. The 
magnetization of the amorphous Ising spins, i.e. the glassy metallic magnet is 
investigated by our modeling and calculation on a theoretical base. The results 
show that the glassy metallic magnets has a lower Curie temperature, weaker 
magnetization, and higher entropy comparing to the regular ferromagnet in 
crystal form. These findings can be understood with the randomness of the 
amorphous system, and agrees well with others’ experimental observations. 

Keywords: Ferromagnetic Metallic Glass, Spontaneous Magnetization, 

Husimi Rhombus Lattice, Structural Glass 


1. Introduction 

Since the first report of Fe-based bulk metallic glass (BMG) in 1995 [H, 
Fe-based BMG with high glass forming ability (GFA) has drawn an intensive 
research interest in the last two decades. The nature of glassy state provides 
advantages of high fracture strength, corrosion resistance, and excellent soft- 
magnetic properties More specifically, difi'erentiated from the conven- 
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tional ferromagnet in the crystal form, the feature of easy magnetization and 
demagnetization of the Fe-based BMG enables various significant applications 
in either scientific research or industrial engineering @ . 

Consequently, in almost every research on BMGs, the study of magnetic 
properties of such materials is an indispensable concern. Some phenomenologi¬ 
cal principles can be indicated with a general review of experimental works, for 
example, the magnetic moments per transition metal atom, the Curie tempera¬ 
ture, and the magnetization of metallic glass are lower than in the correspond¬ 
ing crystal alloy |7Hl2l|: The smaller GFA corresponds to the higher transition 
temperatures, especially the melting transition Tm, i.e. a more thermally stable 
system 12l4l4l| ; And external magnetic field can enhance the soft magnetic prop¬ 


erties of Fe-based bulk amorphous materials [l^ • These observations accord well 
to our intuitive knowledge of the amorphousness. Nevertheless, investigations 
on the featured magnetism of BMG have rarely been conducted on a theoret¬ 
ical base, especially from the aspect of thermodynamics with the fundamental 
concepts of the glass and crystal. 

On the other hand, along with the experimental work, the simulation and 
computation on the metallic glass have also became an activate field [AES- 
However, by this time the theoretical analysis mainly focus on the glass forma¬ 
tion and simple thermodynamics, simulation work involves the concept of spin 
to study the magnetic property is rare. Ref. [l^ as an example. 

According to the above concern, in this work we conduct a theoretical ex¬ 
ploration on the magnetic properties of ferro-BMG with an abstractive thermo¬ 
dynamic model of spins’ crystal and glass. The classical Ising spins are modeled 
on a random-angled recursive lattice to present the magnetic moments arranged 
in an amorphous configuration. The integration of Ising spins and the struc¬ 
tural glass model well represents the case of ferromagnetic metallic glass, while 
the recursive lattice technique enables exact thermodynamic calculation of the 
model. In this way, the spontaneous magnetization transition with Gurie tem¬ 
perature Tc, the effects of randomness parameter, and the entropy behavior are 
quantitatively investigated. 


2. The Model 


Ising model was originally designed to describe the magnetization, and has 

30 . 


been developed to study the magnetic materials for many years 


The 


model is usually applied on lattices for numerical treatment. It should be ad¬ 
dressed that, the model in this paper is not the spin glass 2l|: the orientation 
of spins only relates to the degree of magnetization, and the amorphousness 
of glass is represented by the configurational randomness, i.e. the system is a 
structural glass with magnetic spins as particles. 


2.1. Lattice Construction 

Among numerous methods developed to solve the Ising model, the recursive 
lattice has been proved to be a reliable method with the advantages of simple 
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Figure 1: The recursive lattices constructed from regular lattices, fixed (A & 
A’) and random (B & B’) conformational respectively. The fixed recursive one 
is the classical Husimi lattice (A’). The lattice B’ used in this paper employs 
rhombus units with variable angles to mimic the random lattice B. Note that 
in A’ and B’ all the edges are in the same length, however it cannot be drawn 
on a page in that way. Four spins are presented in each lattice to demonstrate 
the uniform magnetic moment orientation regardless of the structural ordering. 


formulation and exact calculation 22-2^. Recursive lattices is constructed by 
connecting the units from regular lattice on the vertex in a recursive fashion, 
to approximate the regular lattice with the same coordination number. Such a 
lattice assembled by square units is called Husimi lattice [l^. Fig (T] illustrates 
how a Husimi lattice (A’) is drawn from a regular square lattice (A). 

For the structural glass (Fig[Tj3), a calculable regular lattice is impossible. 
However, the unit cell in a recursive lattice only interacts to others on the joint 
vertex and its conformation is independent, this makes a individually random 
unit with variable angles possible. Herein, we construct a Husimi rhombus 
lattice with random angles to demonstrate the conformational stochasticity of 
structural glass (Fig[Tj3’). In this random lattice, the conformation of each rhom¬ 
bus is randomized and fixed by the program, therefore locally it has a quenched 
randomness. 

By setting the square to be the lowest energy conformation, any randomized 
off-90° rhombus will have a higher energy, which implies a more disordered 
conformation deviated from the crystal. Depending on the setup of parameter, 
we can artificially control how “far away” the structure is from the “pure crystal” 
form. Thereafter, it is better to define our amorphous model as “semi-glass”, 
since it is a metastable state other than the crystal, but not reaching the upper 
limit of amorphousness, i.e. the “pure glass”. 
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2.2. Ising model and calculation 

The Ising spins are applied on the lattices as shown in Fig[TJ The structural 
randomness does not affect the arrangement and orientation of the magnetic 
spins, therefore the metallic glassy magnets is possible. In a rhombus unit we 
consider the interactions J between two neighbor spins and Jp between the 
diagonal spin pairs. The Hamiltonian of the system is given by 

e = -jY, -JpY. ( 1 ) 

<i,j> 


where < Si,Sj > is the neighbor pair, < Sj,Sj' > is the diagonal pair, H is 
magnetic field, 9i is the angle against the diagonal pair, and A is the angle 
energy coefficient. 

Assume that J is fixed and the angle only affects the diagonal interaction, 
then by adjusting Jp we can have a system with particular randomness away 
from the ideal crystal state. And further assume that the deviation only affect 
the diagonal interaction in the rhombus due to the diagonal length change, 
and it is linearly proportional to the diagonal length, the we have the diagonal 
interaction energy: 


Ep{e) = -jp 


SjSg 


\/2(r^^cos^ 

\/ 2 (l - cosdo) ’ 


( 2 ) 


where do is the ground state angle with the lowest energy, it is set to be 7r/2 in 
this paper. Then we have A{9) as the angle energy coefficient: 


A[9) = 


1 — cos 9 
1 — cos 9o ’ 


(3) 


Imagine the entire lattice has an original site Sq on the 0th level, then the 
sites nearest to are marked as Si (or S'l to distinguish the two sites in the 
same unit), the site diagonal to Sq is marked as S 2 , and so on. On one rhombus 
unit there are three sub-trees Tn+i and Tn +2 contribnting to the sites Sn, 5"„+i 
and Sn+ 2 , and a larger sub-tree Tn is synthesized by adding the weights of the 
local unit. With the particular spin state of Sn = ±1, define the partial partition 
function (PPF) Z„(S'„) on level n to count the contribution of sub-tree to the 
site Sn (not included), then the nth level PPF is the function of PPFs on the 
previous level n -I- 1, n -I- 2 and the weights of local unit ^( 7 ): 


Zn{ + ) — '^^.Zn+l{Sn+l)Zn+l{Sn+i)Zn+ 2 {Sn+ 2 )w{'^), 


(4) 


7=1 

16 


Zn{ ) — ^ ] Zn+l{Sn+l)Zn+l{Sn^.i)Zn+2{Sn+2)wi'j). (5) 

7—9 

And the partition function (PF) of the whole sub-tree on nth level is given by 

= Zn(+)e^^^- + Zni-)e-^^^\ (6) 
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which summers all the weights at level n including the spin Sn- 

Then this model can be exactly solved by the recursive relations between 
the PF and PPFs on successive levels to obtain the probability of spins’ state 
at particular temperature, and the magnetization and thermodynamics can be 
derived from the probability. For the calculation details please refer to previous 
works on the Ising model on Husimi lattice 2^ 2^3^ and the supplemental 
materials. 


3. Results and Discussion 

In experimental the magnetization of BMG is usuall y in duced by an exter¬ 
nal magnetic field, nevertheless in the 2D Ising model |3l| . spins can sponta¬ 
neously transit to uniform orientation without external field. And the temper¬ 
ature where this spontaneous magnetization occurs can be taken as the Curie 
temperature [^. In this section, all the discussions on the magnetization and 
thermodynamics refer to the case of spontaneous magnetization. 

3.1. The magnetization and Tc 

While the magnetization is usually characterized by the magnetic field strength 
in the real materials, in the theoretical model of spins ensemble the magneti¬ 
zation can be simply presented by the probability of spins orientation. This 
probability of Ising spins can be solved as the ‘solution’ of the model. The 
solutions on fixed Husimi lattice (crystal) and random lattice (semi-glass) with 
J = 1 and Jp = 0.2 are shown in the Fig. [2] Note that the temperature in this 
work is normalized with Boltzmann constant kp = 1- Several observations can 
be drawn here: 1) At high temperature, one spin has fifty percent probability 
to be in either -|-1 or —1 state, then the average orientation is 0; 2) With the 
cooling process, there occurs the spontaneous magnetization in both systems, 
and the Curie temperature is Tq = 3.131 and 3.097 for crystal and semi-glass; 
3) Below Tq, the magnetization continue to increase with the decrease of tem¬ 
perature, and at zero temperature, all the spins will be in the same state and 
we then have a 100% magnetization on the lattice (not shown in the Figd]). It 
is clear that the semi-glass system is less stable and has a lower Tc, i.e. the 
system is easier to proceed a transition, and the magnetization of semi-glass is 
weaker at a particular temperature comparing to the crystal. 

On the other hand, the Tc reduction is relatively small comparing the ex¬ 
perimental results. This is because that our model represents a semi-glass as 
mentioned above, and the parameter Jp = 0.2 in Fig|2] is merely a reference 
setup. We can arbitrarily obtain a particular reduction to fit the real system by 
adjusting the parameter, the details will be discussed below. 

3.2. Paramagnetic and ferromagnetic behavior with external field H 

In experimental, the magnetization responding to the applied external field is 
a critical characterization, e.g. the VSM test, to reveal the magnetic properties. 
In our theoretical model, the field H in EqUcan serve as an adjustment option 
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Figure 2: The degree of magnetization M of Ising spins on random-angled 
Husimi lattice comparing to the regular Ising magnets on fixed lattice, with 
J = 1 and Jp = 0.2. The temperature is normalized with Boltzmann constant 
fcs = 1. 


for the same purpose. Figure [2^ presents the 3D mapping of the magnetization 
with the variation of external field H and temperature T. The paramagnetic and 
ferromagnetic ranges differentiated by the critical temperature can be observed, 
the sign of magnetization indicates the orientation of spins magnetism, and the 
bold discontinuous line is the spontaneous magnetization with H = 0 shown in 
FigH 

Three magnetization curves at T = 3.5, 2.9 and 2.0 are extracted out and 
presented in FiglSb, c and d. While in FiglJJj the paramagnetic behavior is clear 
at high T as expected, the magnetic behavior below Tc is somehow bizarre, that 
the spontaneous magnetization is higher than the system with small field ap¬ 
plied [H = 0.5) in FigjHJ:, and at lower temperature this phenomenon is more 
distinct as shown in FiglSji. The exact reason of this phenomenon is hypothe¬ 
sized as a defective feature of our semi-glass model: Since the model is static, 
the hysteresis loop cannot be simulated, and for the ferromagnetic case we can 
only infer the coercivity by the free energy calculation of self-magnetization case 
with H — 0. While in FiglSJ: and d the extrapolations with H —>■ 0 have a clear 
trend pointing to the origin, implying a soft magnetic behavior even below Tc, 
which is how a real magnetic metallic glass exhibits. Nevertheless, nor like the 
real system where the conformational randomness can cancel out the orientation 
of magnetic spins, the off-crystal structure of the 2D lattice is fundamentally in¬ 
dependent of the spins orientation, subsequently the amorphousness presented 
by random conformation cannot cancel the magnetization. In another word, 
regardless of how the rhombus unit is twisted with higher angling energy, the 
strong coupling between spins at low T is not affected. Therefore, at low tem¬ 
perature the intrinsic magnetism from spontaneously oriented spins will resist 
the soft magnetic behavior and retains a firm coercivity, while the application 
of external field (even a small value) will dominate the spins orientation and 
simulate the real behavior of metallic glass in a field. 
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Figure 3: (a) The 3D mapping of magnetizations with various external field H 
and temperature T. (b, c, d) The magnetization curve at T = 3.5, 2.9 and 2.0. 
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Figure 4: The variation of magnetization M with the ^ 0 in the near-zero 
range in the case of T = 2.9 (Fig.3c). The system shows exact soft magnetism 
with the presence of H, while the self-magnetization case gives a meta-stable 
solutions with H = 0. 


To illustrate the soft magnetism with the presence of the external magnetic 
field, a detailed M vs. H behavior of the near-zero range in the case of FiglSj: 
is presented in the FigHl It is clear to observe that the i? —>■ 0 corresponds 
to M —>■ 0, as an extreme value oi H — 1 x 10“® raises M = 9.84 x 10“^. 
However the calculation of pure self-magnetization case distinguishably gives 
the solutions out of the extrapolation, which is the M = ±34.2% with H — 0 m 
FiglU Notwithstanding, besides this defect, the application of external field on 
this model to simulate the magnetization test is still practical to describe the 
magnetic metallic glass. 

3.3. The thermodynamics 

While various electromagnetic and thermodynamic properties of the glassy 
ferromagnet are being actively researched, the study and report on entropy 
is rare. A probable reason of this lack of attention maybe due to that the 
entropy is not a directly measurable variable. Nevertheless, considering that 
this research field is fundamentally based on the natural difference between 
crystal and amorphous state, should entropy be an important quantity to pay 
attention on. Some experimental work reported the positive correlation between 
higher entropy with enhanced GFA or magnetization . 

Unlike experimental, in the theoretical calculation the entropy is easy to be 
monitored. Figure[5K shows the comparison between the entropy of Ising magnet 
on random rhombus lattice and in the crystal form around the Tq region. For 
both cases, the transitions on entropy at Tq are clear, and accords well to their 
spontaneous magnetization revealed in Figl^l 

In the Hamiltonian (Eq[T]), the diagonal interaction parameter Jp is the only 
variable related to the randomness of angle, and Jp also has positive value setup 
in the model like J, that is, higher Jp prefers same spins orientation and cancels 
the effect of randomness, therefore we can take Jp to be the inverse randomness 
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Figure 5: (a) The entropy of Ising magnet on the random-angled Husimi lat¬ 
tice (semi-glass) comparing to the regular Ising magnet (crystal), (b) Entropy 
variation with diagonal interaction parameter Jp. 


parameter, i.e. the lower the Jp is, the system is more amorphous and farther 
away from the crystal. The effects of Jp variation on the Tq and entropy are 
shown in Fig. [^5. While all the systems exhibit similar behaviors, the variation 
of Jp by 0.1 from 0.1 to 0.5 give the Tc as 2.990, 3.097, 3.296, 3.401, and 3.598 
respectively, exhibiting a rough linear relationship. It is easy to observe that at a 
particular temperature, the increase of Jp corresponds to a more stable system 
with lower GFA, smaller entropy and higher magnetization. This theoretical 
finding agrees well with others’ experimental observations. 

4. Conclusion 

The Ising spins were modeled on a Husimi rhombus lattice with random 
angles to describe the ferromagnetic metallic glass. A diagonal interaction pa¬ 
rameter Jp related to the angle randomness was added into the Hamiltonian 
to disturb the fixed configuration, subsequently provided a semi-glass deviated 
from ideal crystal. The model was exactly solved and exhibits the spontaneous 
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magnetization with transition temperature Tq- The glassy system presents a 
lower Tq and weaker magnetization than the ferromagnetic crystal. Both para¬ 
magnetic and ferromagnetic behaviors at various T with external field to mimic 
the magnetization test were investigated, with a general expectable observa¬ 
tions, that the system behaves as ideal soft-magnetism with the presence of 
external field H, however a higher spontaneous magnetization was also found as 
a meta-stable solution of the pure self-magnetization system with H = 0. The 
entropy and the effect of parameter Jp were studied, the observation agrees 
well with others’ experimental findings, that the thermodynamically less stable 
system corresponds to a higher GFA, larger entropy and weaker magnetization. 
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